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We numerically demonstrate that collective bifurcations in two-dimensional lattices of 
locally coupled logistic maps share most of the defining features of equilibrium second-order 
phase transitions. Our simulations suggest that these transitions between distinct collective 
dynamical regimes belong to the universality class of Miller and Huse model with synchronous 
update [Marcq et al., Phys. Rev. Lett. 77 (1996), 4003]. 

§1. Introduction 

For strong enough coupling, mutually interacting oscillators tend to synchronize. 
Prof. Kuramoto first gave a solid mathematical grounding to this rather intuitive 
idea in 1975. 1 ) More precisely, he showed analytically that an ensemble of non- 
identical phase oscillators with distributed natural frequencies, and coupled through 
their mean-field, synchronize above a critical value of the coupling constant. 2 ) This 
dynamical phase transition is a robust phenomenon, generically occurring in a wide 
class of similar coupled dynamical systems. 3 ) 

The system of interest in the present work is complementary to Kuramoto's 
model: N identical units with aperiodic individual (discrete time) dynamics sit on 
the nodes of a regular lattice, and are updated synchronously with local (diffusive) 
coupling. 4 ^ For a large enough coupling constant, the spatially-averaged activity of 
such coupled map lattices is time-periodic. 5 ) However, synchronization of individual 
units is not involved, at least in the usual sense: macroscopic coherence coexists with 
microscopic disorder, as evidenced by broad distribution functions of local activity 
and rapidly decaying spatial correlation functions. Extensive numerical simulations 
have confirmed that these emergent macroscopic cycles are global attractors, well- 
defined in the infinite-time, infinite-size (thermodynamic) limit: 5 ) fluctuations of the 
spatially-averaged activity about the collective cycle vanish when N — > 00. 

For a large enough control parameter, a single logistic map exhibits an inverse 
bifurcation cascade between regimes of banded chaos. In the same parameter region, 
and for strong enough coupling, locally coupled logistic maps exhibit an inverse 
bifurcation cascade between macroscopic cycles of period 2 ra . 5 ) Generalized mean- 
field arguments give a satisfactory understanding of the build-up of correlations 
at the origin of the dynamical long-range order involved here, 6 ) at least far from 
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bifurcation points. 

In this work, we wish to characterize, thanks to numerical simulations, the first 
bifurcation points in the cascade. 7 ) Upon defining adequate order parameters, we 
demonstrate that these macroscopic bifurcations harbor the characteristic features 
of equilibrium second-order phase transitions. The relevance of standard finite-size 
scaling theory indicates that they persist in the infinite-size limit. Many properties 
of an equilibrium system close to a second-order transition turn out to be largely 
independent of the microscopic details of the interactions between individual compo- 
nents: they fall instead into a small number of universality classes, each defined by 
global features such as the symmetries of the underlying Hamiltonian or the spatial 
dimensionnality of the system. We numerically evaluate the critical exponents of 
period-doubling phase transitions, and discuss the relevance of the notion of univer- 
sality to temporal spontaneous symmetry breaking. 

§2. Critical properties of a period-doubling phase transition 

We consider the dynamics of a set of N = 1? variables x\ j sitting on the nodes 
of a two-dimensional square lattice. Time is discrete. The update rule consists of 
two stages: each site is first updated according to the logistic map with parameter 
r € [0,4]: 

/: [0,1] -[0,1] 

x i— > r x (1 — x), 

then transformed according to a diffusive coupling operator. Interaction is restricted 
to nearest-neighbors. All sites are updated synchronously according to the rule: 

xg 1 = (1 - Ag) f{x\J) + g + /«,_!) + f(xl +1J ) + , (2-2) 

where g is the coupling constant, set to g = 0.2 in this Section (democratic coupling). 
The microscopic control parameter of the lattice dynamical system is r. 
The bifurcation diagram of the mean activity 

1 L 

xt = 12 E x h ( 2 -3) 

is given in Fig. ^ for r > = 3.57. . .. The statistical behavior of the coupled 
map lattice may be interpreted as long-range order accompanied by the temporal 
evolution of spatially-averaged quantities: fixed point (or period 1) above ri_2 ~ 
3.86, period 2 for r € [?"2_4, with r2_4 — 3.63, period 2 n , n > 2, for r S 

[f"oo> r 2-4] • An infinite cascade of period-doubling bifurcations is expected to occur 
in the limit r — > r^. 8 ) 

For simplicity, we focus on the period 1-period 2 bifurcation. The (time-asymptotic) 
distribution of site values in these two regimes is shown for typical parameters in 
Fig. [2j The Lyapunov spectrum scales linearly with the system size: chaos is ex- 
tensive. Dynamical quantifiers such as the largest Lyapunov exponent and the Kol- 
mogorov entropy remain continuous close to the bifurcation point (see also 9 )). 
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Fig. 1. Bifurcation diagram of democratically coupled (g = 0.2) logistic maps on a two-dimensional 
lattice of L 2 — 1024 2 sites: for each parameter value r > roo, 20 consecutive values of the mean 
activity x* are plotted vs. r. Initial site values are randomly distributed over the interval [0, 1], 
and a transient of duration to = 10 4 is discarded. Period-doubling bifurcations are rounded by 
finite-time, finite-size effects. 




Fig. 2. Probability distribution function of the lattice's site values, measured for g = 0.2, L = 1024, 
in the period 1 (left-hand graph, r = 3.93) and period 2 (right-hand graph, r — 3.75) collective 
states. The pdfs measured over 4 consecutive time-steps are superposed. The breadth of the 
lines is related to finite-size statistical fluctuations. 



The local "magnetization" is defined by rn\_ 2 = x — x . The order parameter 
of the transition reads: 7 )' 10 ) 

Mi_ 2 = ( K_ 2 | ) = T lim I W t+1 ~ * 2t \ • (2-4) 

The period 1 and period 2 phases are respectively "paramagnetic" (Afi-2 = 0) and 
"ferromagnetic" {M\-2 ^ 0). Through Eq. Q2-4j) . we wish to draw a strong analogy 
with the Ising model. Instead of spin-reversal invariance, time-translation invariance 
is broken in the ordered phase of the lattice dynamical system. Further, the correla- 
tion length £ exhibits a sharp maximum close to the transition, where £ characterizes 
the exponential decay of the equal-time, two-point correlation function of the field 
x\y Large scale simulations suggest the presence of well-defined power laws control- 
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Fig. 3. Period 1-period 2 transition, g — 0.2, r = rjl 2 = 3.86212. The finite-size scaling laws 
given in Eq. I2-5H allow to measure the exponent ratios (5/v (a), j/v (b), \jv (c). Corrections 
to scaling are apparent in the three graphs. The dashed lines correspond to the exact exponent 
values of the two-dimensional Ising model. 
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The precise location of the transition point is determined using Binder's method: 11 ) 



the value of the finite-size cumulant J7fl 2 ( r ) = — 3 + \(rn\_ 2 ) ) / (("1-1-2) / a ^ 
size N = L 2 is independent of L at criticality: VL, Ui_ 2 {r^_ 2 ) = U^ 2 . Us- 
ing system sizes ranging from L = 32 to L = 128, we find r^_ 2 = 3.86212(12), 
—U^t 2 = 1-819(12), where numbers within parentheses correspond to the uncer- 
tainty over the last digit(s). The exponent ratios /3/v, "f/v, 1/f are then obtained 
from the following scaling laws: 



(2-5) 



Mi_ 2 ( J L) 
Xi-2(L) 
d r U[_ 2 

9 r logMi_ 2 (L) 
d T \og{m\__ 2 (L) 2 ) 

where statistical averages are computed at the infinite-size critical parameter value 
r i^2- For all the above observables 0, our data is consistent with the simplest 
equilibrium correction to scaling: O = (ao + aiL~ w ° + ...), where 4>o and 000 
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are respectively the scaling exponent denned in Eqs. (|2 • 5(1 and the first subdominant 
exponent. Taking into account this correction, our best estimates are (see 12 ) for a 
detailed account of the procedure we use) : 

P/u = 0.126(4), 

1 /u = 1.76(3), (2-6) 
v = 0.865(25). 

The exponent ratios (3/v and 7/V are in excellent agreement with the two- 
dimensional equilibrium Ising model: (/3/V)i s ing = 1/8, ^/i^ising = 7/4. However, 
the correlation length exponent u is not consistent with z^i s i ng = 1 (see Fig. |2J) . As a 
consequence, our estimates for (3 = 0.108(6) and 7 = 1.52(7) are also incompatible 
with the Ising values. Unless corrections to scaling of an unusual nature are present 
in this system, our data suggests that the period 1-period 2 phase transition does 
not belong to the universality class of the two-dimensional equilibrium Ising model: 
the analogy we have drawn breaks down at a quantitative level. 

§3. Universality 

The evolution rule (|2-1() - (I2-2|) involves two parameters: the coupling constant g 
and the nonlinear parameter r: a line of period 1-period 2 transitions is observed 
in the parameter plane for large enough coupling (g > 0.10). We studied two other 
transition points on this line: a r-driven transition at g = 0.11, and a coupling-driven 
transition at a fixed value of r = 3.831. The same measurement protocol is used, 
with similar system sizes and statistical accuracy (see 7 ) for further details). In the 
three cases, we find exponent values mutually consistent within error bars (see Table 
for a summary of results). This suggests that the period 1-period 2 transition line 
defines a universality class distinct from that of the two-dimensional Ising model. 
However, the hyperscaling relation 2(5 + 7 = dv [d = 2) remains valid, and the 
critical value U°° of Binder's cumulant is consistent with the Ising value. 

We also studied transitions between cycles of higher period, which are charac- 
terized by the same phenomenology as described above. Physical observables are 
defined as before, upon replacing m'_ 2 by the appropriate local magnetization, for 
instance ra\_^ = x 4 ' +2 — x 4 * for period 2-period 4 transitions. We find that: (i) these 
transitions do not belong to the Ising universality class; (ii) the measured exponent 
values are consistent with those found for period 1-period 2 transitions - even though 
the system sizes we use do not allow a clear-cut answer to this last question. 

In fact, the same features are characteristic of other ordering transitions between 
chaotic phases. Miller and Huse introduced a two-dimensional coupled map lattice 
with a microscopic "up-down" symmetry, using a continuous, piecewise-linear, odd- 
symmetric map. 13 - 1 In a previous work, 12 ) we showed that the Ising-like transition 
of this extensively chaotic lattice dynamical system does not belong to the Ising 
universality class, with in particular v = 0.89 ± 0.03. However, Ising exponents are 
recovered for the same geometry, coupling, and local map once the update rule be- 
comes asynchronous (lattice sites updated one at a time). The nature of update, a 
dynamical feature, is a relevant "parameter" that distinguishes between universality 
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2D Ising 


MH 


5 = 0.2 


5 = 0.11 


r = 3.831 


Threshold 




0.20534(2) 


3.86212(12) 


3.8310(1) 


0.1100(1) 


-U°° 


1.83 


1.832(4) 


1.819(12) 


1.828(10) 


1.83(3) 


PI" 


0.125 


0.125(4) 


0.126(4) 


0.131(10) 


0.129(10) 


l/v 


1.75 


1.748(10) 


1.76(3) 


1.75(4) 


1.73(3) 


(2/3 + i)/u 


2 


2.00(2) 


2.02(4) 


2.01(6) 


1.99(5) 


P 


0.125 


0.111(5) 


0.108(6) 


0.107(21) 


0.107(10) 


7 


1.75 


1.55(4) 


1.52(7) 


1.50(12) 


1.39(10) 


V 


1 


0.887(18) 


0.865(25) 


0.860(45) 


0.80(4) 



Table I. Summary of critical exponents: Ising model; Miller and Huse model; r-driven period 1- 
period 2 transition at g — 0.2; r-driven period 1-period 2 transition at g — 0.11; g-driven period 
1-period 2 transition at r — 3.831. 

classes. The critical exponents of Miller and Huse's model with synchronous update 
are recalled in the Table: we showed 12 ) that they are representative of a genuine 
universality class in the sense that changes of, e.g., the local map do not alter them 
provided that the up-down symmetry is preserved (see also 14 -*' 15 -*). Remarkably, 
these exponents are consistent (within error bars) with those obtained for the period 
1-period 2 transitions of two-dimensional coupled logistic maps. The same universal- 
ity class encompasses ordering transitions between chaotic phases of synchronously- 
updated lattice dynamical systems, whether due to the (discrete) breaking of an 
up-down symmetry in phase space or to that of time-translation invariance. 

§4. Conclusion 

Our numerical simulations demonstrate that equilibrium finite-size scaling laws 
allow to characterize the period-doubling macroscopic bifurcations of lattices of lo- 
cally coupled chaotic logistic maps. This relevance provides in itself further evidence 
that such bifurcations are indeed continuous phase transitions, well-defined in the 
infinite-size limit, where the equal-time correlation length diverges as the system's 
linear size. We evaluate the static critical exponents /3, 7 and v of fixed point-period 
2 and period 2-period 4 transitions. The best interpretation of our (finite-size) data 
is the following: these transitions belong to the universality class of the Miller and 
Huse model with synchronous update, not to the equilibrium Ising universality class: 

V ~ 0.89 / Rising = 1. 

Period-doubling phase transitions are by all means unusual: they separate mi- 
croscopically chaotic states with different (regular) macroscopic dynamics. The ex- 
istence of an ordered phase is related to the discrete breaking of time-translation 
invariance at the macro-scale, while microscopic dynamical quantifiers, such as the 
Lyapunov spectrum, remain continuous close to the transition. The lattice dynamical 
system is a priori far from equilibrium: microreversibility is not expected to hold. 
The non-Ising value of the correlation length exponent v, at odds with standard 
coarse-graining arguments, 13 )' 16 ) remains a puzzle. 

Two comments are in order. First, the control parameters we use are micro- 
scopic. This is also true of related models with similar critical properties. 14 )' 15 )' 17 ) 
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The various phases we observe arise due to a balance between nearest-neighbor in- 
teractions and fluctuations of deterministic origin. Since the critical exponents we 
measure are not those of mean-field theory, 18 ) we know that fluctuations cannot 
be neglected close to the transition points. One would certainly like to be able to 
satisfactorily define the temperature of a given extensively chaotic lattice dynami- 
cal system, i.e. a macroscopic, intensive parameter, well-defined in the infinite-size 
limit, that quantifies at a coarse-grained level the degree of microscopic disorder. 
This remains an open, challenging problem. 

Second, the synchronous nature of the update rule of a coupled map lattice is 
never insignificant. 19 )> 20 ) In the case of Miller and Huse model, Ising exponents are 
recovered with asynchronous update: the nature of update is a relevant "parameter" 
in the sense of critical phenomena. For coupled logistic maps, the long-range order 
leading to periodic collective behavior is destroyed by asynchronous update. The 
fully synchronized state, where all lattice site values are equal to that of the unstable 
fixed point of the local map, remains unstable under synchronous update. However, 
synchronization transitions occur as soon as the update rule includes some degree of 
asynchrony. 21 ) In all cases, the choice of a given update rule is a crucial modeling 
issue. 
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